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Abstract 

Background: The atrial gap junction protein connexin-40 (Cx40) has been implicated to play an important role in 
atrial conduction and development of atrial fibrillation (AF). However, the frequency of Cx40 mutations in AF 
populations and their impact on Cx40 expression remains unclear. In this study, we sought to identify 
polymorphisms in the Cx40 gene GJA5, investigate the potential functional role of these polymorphisms, and 
determine their allelic frequencies. The prevalence of nonsynonymous Cx40 mutations in blood and atrial tissue 
was also compared to mutation frequencies reported in prior studies. 

Methods: We conducted direct sequencing of the GJA5 coding and 3' UTR regions in blood samples from 91 lone 
AF subjects and 67 atrial tissue-derived samples from a lone cohort, a mixed AF cohort, and several transplant 
donors. Reporter gene transfection and tissue allelic expression imbalance assays were used to assess the effects of 
a common insertion/deletion polymorphism on Cx40 mRNA stability and expression. 

Results: We identified one novel synonymous SNP in blood-derived DNA from a lone AF subject. In atrial 
tissue-derived DNA from lone and mixed AF subjects, we observed one novel nonsynonymous SNP, one rare 
previously reported synonymous SNP, and one novel 3' UTR SNP. A previously noted 25 bp insertion/deletion 
polymorphism in the 3' UTR was found to be common (minor allele frequency = 0.45) but had no effect on Cx40 
mRNA stability and expression. The observed prevalence of nonsynonymous Cx40 mutations in atrial tissues derived 
from lone AF subjects differed significantly (p = 0.03) from a prior atrial tissue study reporting a high mutation 
frequency in a group of highly selected young lone AF subjects. 

Conclusions: Our results suggest that Cx40 coding SNPs are uncommon in AF populations, although rare 
mutations in this gene may certainly lead to AF pathogenesis. Furthermore, a common insertion/deletion 
polymorphism in the Cx40 3' UTR does not appear to play a role in modulating Cx40 mRNA levels. 
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Background 

Atrial fibrillation (AF) is the most commonly encoun- 
tered form of sustained cardiac arrhythmia in clinical 
practice, with a prevalence in the USA that is projected 
to increase three-fold between the years 2000 and 2050 
[1-3]. AF is a complex disease characterized by irregular 
electrical activity within the atria, resulting in episodes 
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of uncoordinated atrial contraction that substantially in- 
crease stroke risk and mortality [4-6]. Studies investigat- 
ing the genetic basis of the disease have shown that 
somatic and germ line variants in the atrial gap junction 
protein connexin-40 (Cx40) can disrupt normal atrial 
conduction and may predispose individuals to idiopathic, 
or lone, AF [7-9]. The Cx40 gene GJAS encodes two al- 
ternative transcripts, Cx40-A and Cx40-B, comprised of 
first exons 1A or IB and a shared coding second exon, 
respectively [10]. We previously described rs 10465885, a 
common single nucleotide polymorphism (SNP) in the 
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Cx40-B promoter that strongly affects Cx40 mRNA ex- 
pression and is associated with early onset lone AF [11]. 
We hypothesized that direct sequencing of the GJAS re- 
gion in atrial and blood DNA from expanded cohorts of 
mixed and lone AF patients might yield additional poly- 
morphisms altering Cx40 structure or expression, poten- 
tially resulting in a predisposition to AF. 

Methods 

Study subjects 

Most (64 of 67) of the atrial appendage tissues were 
obtained from a biorepository of human atrial tissues from 
patients undergoing coronary artery bypass grafting, valve 
surgery and/or the Maze surgical ablation procedure. Sub- 
jects had a history of AF and were categorized as "lone 
AF" if they did not have structural heart disease, or "mixed 
AF" if they had AF in combination with coronary artery 
disease or valvular heart disease (mitral regurgitation). 
The remaining three atrial tissue samples were obtained 
from non-failing donor hearts not used for transplant with 
unknown cardiac phenotypes (Table 1). Use of discarded 
atrial tissues from surgical subjects was approved with 
an exemption for written consent by the Cleveland 
Clinic Institutional Review Board and performed in ac- 
cordance with institutional guidelines. Atrial specimens 
were obtained from unmatched organ donor hearts, 
after informed consent by a surviving family member 



and used under a protocol approved by the Cleveland 
Clinic Institutional Review Board. 

Blood samples were obtained from the Cleveland 
Clinic Lone Atrial Fibrillation GeneBank (CCAF), which 
enrolled subjects at least 18 years of age who had a his- 
tory of AF without coronary heart disease. Specifically, 
these subjects had no significant coronary artery disease 
(<50% coronary artery stenosis if cardiac catheterization 
was performed or a normal stress test) and normal left 
ventricular function (left ventricular ejection fraction 
>50%). Patients were excluded if they had structural 
heart disease or congenital heart disease, except for an 
isolated patent foramen ovale with normal right heart 
chambers; significant valvular disease (>2+ regurgitation 
or any valve stenosis); or prior percutaneous coronary 
intervention or coronary artery bypass grafting. Patients 
with hypertension were not excluded. Written informed 
consent was obtained from all CCAF subjects under a 
protocol approved by the Cleveland Clinic Institutional 
Review Board and performed in accordance with institu- 
tional guidelines. There were no subjects from whom we 
obtained both atrial tissue- and blood- derived DNA. 

Preparation of genomic DNA, total RNA and 
complementary DNA from human atrial tissue and blood 

Atrial genomic DNA (gDNA) and total RNA was pre- 
pared from atrial tissue samples using Qiagen kits. 



Table 1 Patient characteristics 







Atrial tissues 




Blood samples 




Lone AF 


Mixed AF 


Donor 


lone AF 




N=34 (%) 


N=30 (%) 


N=3 (%) 


N=91 (%) 


Age, years, mean ± SD (range) 


53.8 ± 9.6 


66.3 ± 9.8 


50 ± 15.6 


57.8 ±11.5 




(31-71) 


(32-86) 


(32-60) 


(25-84) 


Sex, male 


28 (82.4) 


1 7 (56.7) 


2(67) 


67 (73.6) 


Race, Caucasian 


32 (94.1) 


29 (96.7) 


3 (100) 


87 (95.6) 


Unknown or Other 


2 (5.9) 


1 (3.3) 


0 


4 (4.4) 


Cardiac Disease 






ND 




CAD 


0 


3 (9.1) 




0 


Valvular disease 


0 


21 (63.6) 




0 


CAD + Valve disease 


0 


6 (18.2) 




0 


Unknown 


0 


0 




0 


AF Burden 






ND 




Paroxysmal 


12 (35.3) 


6 (20.0) 




39 (42.9) 


Persistent 


6 (17.6) 


5 (16.7) 




33 (36.3) 


Permanent 


16 (47.1) 


19 (63.3) 




13 (14.3) 


Unknown 


0 


0 




6 (6.6) 


Rhythm at the time of sampling 






ND 


ND 


Sinus 


16 (47.1) 


10 (33.3) 






AF 


18 (52.9) 


20 (66.7) 






Unknown 


0 


0 







AF = atrial fibrillation, CAD = coronary artery disease, ND = not determined. 



Tchou et al. BMC Medical Genetics 2012, 13:102 
http://www.biomedcentral.com/1471-2350/13/102 



Page 3 of 9 



cDNA was produced with the iScript Select cDNA 
Synthesis Kit (BioRad) using an oligo(dT) 2 o primer. 
Buffy coat gDNA was prepared from blood samples of 
subjects in the CCAF Genebank using the MasterPure 
DNA Purification Kit for Blood Version II (Epicentre 
Biotechnologies). 

Firefly luciferase-G-/>45 3' UTR Expression Vector 
Construction 

The GJAS 3' UTR region (1,938 bp) from atrial tissue 
gDNA was amplified by PCR using the upstream primer 
ATGATCTCGAGAAGCGACGTCTTAGTAAGGCCAG 
and the downstream primer ATGCATGGGCCCCCTT- 
TACCCATCCCATCAGCACC (Xhol and Apal sites in 
bold). PCR products were cloned into the pGEMT-easy 
vector (Promega), transformed into DH5a Escherichia 
coli, and screened for the presence of the insert. Wild- 
type and variant GJAS 3 r UTR insert sequences were 
then directionally cloned into a modified pcDNA3-Luc 
vector, immediately downstream of the luciferase coding 
region. The resulting constructs were co-transfected into 
HL-1 murine atrial cardiomyocytes [12] with a p- 
galactosidase-containing transfection control vector, and 
the Dual-Light reporter gene assay system (Tropix) was 
used to measure luciferase and (3-galactosidase activities 
in cell lysates according to methods previously described 
[11]. The mRNA instability positive control was gener- 
ated by annealing the single-stranded oligos TCGA- 
CAT TT AT TTATAT TT ATT TACCGC GC GGC GCCG 
and T CG ACG G CG CCG CG CG GTA AATA A ATATA AA- 
TAAATG, which contain an AU-rich element (ARE, bold 
emphases) that reduces mRNA stability when present in 
the 3'UTR [13]. Annealed oligos were phosphorylated with 
T4 polynucleotide kinase (M0201, NEB) and cloned into 
the Xhol site of the GJAS 3 r UTR luciferase reporter 
described above. 

GJAS gDNA and cDNA Sequencing 

Screening of 67 atrial derived gDNA samples for the 25 
bp insertion/deletion polymorphism was carried out 
using the PCR primers ATTCCTCGGAGTAGTGGT- 
GAGATGG (upstream) and ACACCCTAGCAGAAG- 
GAAAGGTTGC (downstream), which amplify a 131 bp 
region in the GJAS 3 r UTR. For sequencing of the GJAS 
coding region in gDNA derived from blood or atrial tis- 
sue, two overlapping fragments of GJAS were amplified 
using the primer pairs ACGAGTACCCGGTGGCAGA- 
GAAGGC (upstream) and AGGAGCCAAGCAGTGAT- 
GACAGTGAGAA (downstream) and GCTAATATGG- 
CTACTTTGAATCTTCTC (upstream) and ACATGC- 
AGGGTGGTCAGGAA (downstream), respectively. The 
GJAS PCR fragments were sequenced and analyzed for 
SNPs using the SNP Analysis tool hosted by the La- 
boratory of Population Genetics, National Cancer 



Institute, NIH (http://lpg.nci.nih.gov/LPG) [14]. To 
measure allelic expression imbalance, we used a highly 
sensitive sequencing technique that compares cDNA 
and gDNA allelic ratios in heterozygous individuals 
[15]. This is an adaptation of a similar technique using 
single base extension that also compares allelic ratios in 
cDNA vs. gDNA in heterozygous subjects [16,17]. A 
Cx40-A fragment containing exon 1A, the coding 
region and approximately 1 kb of the 3 r UTR (including 
the 25bp insertion/deletion and an indicator SNP 
rsl043806) was amplified from eight atrial tissue- 
derived cDNA samples using the PCR primers GGTGG 
AAGAGGAACAACTGA (upstream) and GGGCCTCC 
ATAGCTGTCATCA (downstream). A portion of the 
GJAS 3' UTR containing the indicator SNP rsl043806 
was amplified from gDNA using the PCR primers 
ATTCCTCGGAGTAGTGGTGAGATGG (upstream) and 
GGGCCTCCATAGCTGTCATCA (downstream). The 
Cx40-A cDNA fragment and the gDNA GJAS 3' UTR 
fragment corresponding to each subject were sequenced 
using Sanger methodology on an Applied Biosystems 
3730x1 DNA analyzer. The sequence trace files were 
analyzed using PeakPicker 2 software (publicly available 
at http: 1 1 genomequebecmcgillca I publications! pastinen I) 
[15] to determine the cDNA allelic expression ratio, 
normalized to the gDNA ratio, at the rsl043806 pos- 
ition as previously described [11]. 

Results 

Patient Characteristics 

Table 1 shows the demographics of 91 lone AF blood 
donors and a separate group of 67 subjects who pro- 
vided atrial tissue. Of the 67 atrial appendage 
donors, 34 had lone AF (no significant structural 
heart disease), 30 had AF with coronary artery dis- 
ease or valvular heart disease (mitral regurgitation), 
and 3 were derived from heart transplant donors 
with limited phenotypic information. Approximately 
half (56.7%) of the 67 tissue donors were in AF at 
the time of surgery. 

Single nucleotide polymorphisms in the connexin-40 
coding region 

In order to screen for and assess the frequency of nonsy- 
nonymous Cx40 variants in a large cohort of subjects of 
European ancestry, we conducted direct sequencing of 
the Cx40 coding region in blood- and atrial tissue- 
derived gDNA samples from patients diagnosed with 
lone or mixed AF, as well as a small number of trans- 
plant donors. Detected SNPs were determined to be 
novel if they have not been previously reported in public 
SNP databases, including those available from dbSNP 
[18] and the Thousand Genomes Project [19]. Out of 34 
atrial tissue samples obtained from lone AF patients 
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(Table 1), a single novel nonsynonymous Cx40 coding 
SNP was identified (339G>C) (Figure 1A). The novel 
339G>C SNP results in a substitution of serine for argin- 
ine in the cytoplasmic loop of Cx40 (R113S) at a pos- 
ition that is highly conserved in mammals. Analysis with 
the PolyPhen-2 tool predicts that R113S is likely to be 
functionally and structurally well-tolerated [20]. Since 
blood samples matching our atrial tissue samples were 
not available, we were unable to determine whether the 
339G>C variant had a somatic or germline origin. No 
novel SNPs were observed through direct sequencing of 
33 atrial tissue samples from a mixed AF and donor co- 
hort (Table 1). Sequencing of 91 blood samples from a 
lone AF cohort (Table 1) yielded one novel synonymous 
Cx40 coding SNP (951T>C, Figure IB) and one previ- 
ously reported synonymous coding SNP (rs2232191) in 
separate individuals. A search in dbSNP shows that the 
latter SNP is rare in the HapMap European and African 
populations, but has an allele frequency of ~ 5% in the 
Chinese and Japanese populations [18]. 

Identification of variation in the proximal region of the 3 / 
UTR 

We conducted a search for novel regulatory variants of 
GJAS by sequencing the proximal 1,435 bp of the GJAS 
3 f UTR region in atrial tissue genomic DNA samples 
from 48 patients diagnosed with AF, a subset of the 67 
subjects described in Table 1. We observed a novel poly- 
morphism 22 bp downstream of the stop codon 
(*22G>A) in a single sample. We also observed a 25 bp 
insertion/deletion in multiple subjects (Figure 2A-C). 
The insertion/deletion, described in data uploaded by 



Levy and colleagues in the first reported human diploid 
genome sequence, is located approximately 640 bp 
downstream of the stop codon and approximately 1,300 
bp upstream of the poly-A signal [21]. We developed a 
PCR length polymorphism assay to screen for the inser- 
tion/deletion in the full set of atrial gDNA derived from 
67 tissue donors (Figure 2D), and we determined the 
minor allele (insertion) frequency to be 0.45. We also 
frequently observed the common 3' UTR SNP 
rsl692141, located 53 bp downstream of the GJAS stop 
codon with a minor allele frequency of 0.237 in our 
samples. 

Linkage disequilibrium of the GJA5 3' UTR insertion/ 
deletion 

Haplotype analysis of the GJAS 3' UTR insertion/dele- 
tion was conducted with the atrial tissue samples, exam- 
ining eight SNP genotypes situated within the GJAS 
region. The SNP genotypes tested consisted of the com- 
mon Cx40-B promoter variant rsl0465885, a Cx40-A 
promoter SNP (rsl 1552588) and five other contiguous 
GJAS SNPs from the Illumina Hap550 BeadChip 
[11,22,23]. One of these SNPs, rsl043806 in the GJAS 3' 
UTR, was used as the indicator SNP in an expression 
imbalance assay, as described below. Analysis of linkage 
disequilibrium (LD) relationships using Haploview v.4.1 
(Broad Institute) determined that the insertion/deletion 
polymorphism is in partial LD with rsl0465885, with 
r 2 = 0.81 (Figure 3). We have previously shown that 
rs 10465885 is modestly associated with early onset lone 
AF, exerting a strong cis -effect on atrial expression of 
the Cx40-B transcript, and is also associated with atrial 



113 

REAER/SAKEVR 

ACGGGAGGCCGAGAGqGCCAAAGAGGTCCGG 

339 



R 317 

D IQVRYGQKPE 

catccaggttcgttaJggccagaagcctgag 

Figure 1 Sequence traces of novel Cx40 SNPs. (A) Sequence of 339G>C, a novel nonsynonymous Cx40 SNP identified in a lone AF atrial tissue 
sample. 339G>C results in the cytoplasmic loop substitution R113S. (B) Sequence of 951T>C, a novel synonymous Cx40 SNP identified in a lone 
AF blood sample. 
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Figure 2 Identification of a 25 bp insertion/deletion polymorphism in the GJA5 3' UTR. Sequence traces from human atrial gDNA from (A) 
a deletion homozygote, (B) an insertion homozygote, and (C) a heterozygote are shown. Nucleotide sequence numbering is based on the GJA5 
Refseq entry (NM_181703). (D) PCR assay of atrial gDNA from 15 individuals, identifying homozygotes for the insertion (upper band, 206 bp), 
homozygotes for the deletion (lower band, 181 bp), and heterozygotes (with both bands). Screening of 67 atrial samples determined the minor 
allele frequency of the insertion to be 0.45. 



total Cx40 mRNA levels [11]. We hypothesized therefore 
that the insertion/deletion polymorphism might also 
play a role in regulation of Cx40 expression. 

Functional assay of the GJA5 3' UTR insertion/deletion 

We designed a luciferase reporter assay to investigate 
the potential effects of the common insertion/deletion 
on Cx40 mRNA stability and/or translation. GJAS 3 r 
UTR sequences containing the deletion, the insertion, or 
a combination of the insertion and the 3 r UTR variant 
*22G>A were fused to the 3 r end of the firefly luciferase 
gene in a modified pcDNA3 vector [24]. The resulting 
insertion/deletion constructs were transiently transfected 
into HL-1 atrial murine cardiomyocytes and the cell 
lysates were assayed for luciferase activity, normalized 
for transfection efficiency. Comparable levels of activity 
were observed for cells expressing each construct 
(Figure 4A), suggesting that the insertion/deletion poly- 
morphism is not associated with altered Cx40 tran- 
slational efficiency or mRNA stability. As a positive 
control, a luciferase reporter construct containing the 
GJAS 3'UTR sequence with the insertion polymorphism 
was modified to include an AU-rich element (ARE) that 



has been shown to decrease mRNA stability when 
present in the 3'UTR [13]. HL-1 murine cardiomyocytes 
transfected with the ARE- GJAS 3'UTR construct demon- 
strated a 54% reduction in luciferase activity compared 
to cells transfected with the unmodified reporter 
(Figure 4B, p<0.05). 

Allelic expression imbalance assay for the GJA5 3' UTR 
insertion/deletion 

The potential effect of the insertion/deletion on atrial 
Cx40 mRNA levels was also tested using a previously 
described allelic expression imbalance assay [11,16]. In 
this assay, the relative expression of two allelic Cx40 
transcripts was quantified by direct cDNA sequencing in 
subjects heterozygous for a common indicator SNP con- 
tained within the mRNA. PeakPicker 2 software [15,25] 
was used to calculate the allelic ratios by measuring peak 
heights from cDNA sequence traces, normalized to the 
allelic ratio of the corresponding heterozygous gDNA. 
Of the two alternatively spliced Cx40 transcripts, Cx40- 
A mRNA expression was measured in order to exclude 
the strong allelic imbalance observed for the Cx40-B 
transcript in rsl0465885 heterozygotes. We tested eight 
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Figure 3 Linkage disequilibrium of the G-/A5 3' UTR insertion/ 
deletion. Linkage disequilibrium relationships (r 2 ) between the 
Cx40-B promoter SNP (rs1 0465885), a Cx40-A promoter SNP 
(rs1 1552588), the 3' UTR SNP (rs1 043806) used for the allelic 
expression imbalance assay, the insertion deletion polymorphism 
(Ins/Del), and four other GJA5 SNPs from the lllumina Hap550 
BeadChip. The positions of the SNPs relative to the alternative 
transcripts for Cx40-A and Cx40-B are displayed. The D' between 
rsl 043806 and the insertion/deletion is 1.0 (not shown). 



atrial tissue samples heterozygous for the insertion/dele- 
tion and the indicator SNP rs 1043806. The minor G al- 
lele of rs 1043806 is in partial LD with the insertion/ 
deletion (r 2 = 0.31, Figure 3) and is consistently in phase 
with the deletion allele (D' = 1.0). Six of the eight sam- 
ples demonstrated weak allelic expression imbalance, 
yielding a mean 1.56-fold increase in Cx40 expression 
for the deletion-G allele (Figure 5). The remaining two 
samples showed no change in Cx40 levels or imbalance 
in the opposite direction, respectively (Figure 5), yielding 
an overall mean 1.39-fold increase for the deletion-G al- 
lele. The observed pattern of incomplete expression im- 
balance indicates that the insertion/deletion is not 
directly associated with regulation of Cx40 expression; 
however, the insertion/ deletion polymorphism may be in 
partial LD with a separate unidentified modest-strength 
regulatory variant. 




Del Ins Ins + *22G>A 

3' UTR Type 



B 



7.5 n 




Ins ARE-lns 
3" UTR Type 

Figure 4 The 3' UTR insertion/deletion polymorphism does not 
significantly affect GJA5 mRNA stability or translational 
efficiency. (A) HL-1 murine cardiomyocytes transfected with 
vectors containing the firefly luciferase gene fused at its 3' end to 
the GJA5 3' UTR were measured for luciferase activity, normalized 
to the (3-galactosidase activity of the transfection control vector. 
Shown are the mean activities, ± standard deviations. No 
significant differences in activity were observed for constructs 
carrying the deletion (Del, n=8), the insertion (Ins, n=4), or the 
insertion with a 3' UTR point mutation (Ins + *22G>A, n=4) 
(P = 0.256, one-way ANOVA). (B) A GJA5 3' UTR construct carrying 
the insertion was modified to include a mRNA instability sequence 
upstream of the GJA5 3' UTR. The resulting positive control vector 
(ARE-lns) demonstrated reduced luciferase activity compared to 
the unmodified vector (p<0.05, n=3). 



Statistical analyses of Cx40 missense allele frequencies in 
human blood and atrial tissue 

Fisher s exact probability test was used to compare the 
frequencies of nonsynonymous variant Cx40 alleles 
observed in our studies with those reported by groups 
studying similar sample cohorts. We found no statisti- 
cally significant difference (p = 0.56) between the fre- 
quency of missense Cx40 alleles in our 91 lone AF blood 
samples (variantnon-variant alleles = 0:182, 0%) and the 
frequency observed in 218 lone AF blood samples by 



Yang and colleagues (3:433, 0.69%) [26]. However, we 
observed a significant difference (p = 0.03) when com- 
paring Cx40 missense allele frequencies in our 34 lone 
AF tissues (average age at surgery 54 years), in which we 
detected a single nonsynonymous variant (variant:non- 
variant alleles = 1:67, 1.47%), and the frequency observed 
by Gollob et al. in 15 highly selected and young (average 
age of onset 45 years) lone AF tissue donors, in which 
four novel heterozygous Cx40 missense mutations were 
found (4:26, 13.3%) [9]. 
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Figure 5 Allelic expression imbalance analysis of Cx40-A mRNA. 

Shown are eight atrial tissue samples heterozygous for the GJA5 3' 
UTR insertion/deletion polymorphism and the indicator SNP 
rs1 043806. Each bar represents the allelic expression ratio for an 
individual subject, based on peak heights for cDNA/gDNA at the 
indicator SNP site. Significant allelic mRNA imbalance is found when 
the cDNA allelic ratio falls outside of a 95% confidence interval 
created by the gDNA allelic ratios (dotted lines). The G allele of the 
indicator SNP was always found together with the deletion allele of 
the insertion/deletion polymorphism. Six samples demonstrated 
weak allelic expression imbalance of Cx40 favoring the deletion-G 
allele. One sample was not significantly associated with allelic 
expression imbalance and another sample exhibited Cx40 expression 
imbalance in the opposite direction. 

V ) 



Discussion 

We conducted direct sequencing of the GJAS coding re- 
gion in 91 lone AF blood samples and 67 atrial tissue 
samples (34 lone AF, 30 mixed AF, and 3 donors) from 
individuals of predominantly European ancestry, identi- 
fying a novel nonsynonymous SNP (339G>C), a novel 
synonymous SNP (951T>C), and a previously reported 
synonymous SNP (rs2232191). In the 3' UTR we identi- 
fied a single novel variant (*22G>A). We also deter- 
mined that a 25 bp insertion/deletion polymorphism in 
the GJAS 3' UTR is common and in partial LD with the 
Cx40 transcript B promoter SNP rs 10465885, which is 
associated with early onset lone AF and strongly corre- 
lated with Cx40-B and total Cx40 expression [11]. How- 
ever, reporter gene transfection and Cx40-A allelic 
expression imbalance assays indicated that the insertion/ 
deletion is not directly associated with changes in Cx40 
expression. The modest and inconsistent 1.39-fold effect 
observed for the insertion/ deletion in the 8 subjects 
assayed for Cx40-A allelic expression imbalance was 
much weaker than the consistent 3.3-fold effect that 
rsl0465885 has on allelic expression imbalance of the 
Cx40-B transcript [11]. These allelic expression imbal- 
ance results support the hypothesis that the insertion/ 
deletion is in partial LD with a separate modest 



regulatory variant that is responsible for the partial allelic 
imbalance observed. 

The most recent and largest AF genome-wide associ- 
ation study has not identified common variants in or 
near the GJAS gene to be associated with AF at the con- 
servative genome-wide significance threshold [27]. How- 
ever, we have previously identified a common GJAS 
promoter variant that both alters Cx40 expression levels 
and is weakly associated with AF [11]. Since not all of 
the heritability for common traits can be attributed to 
common variants, sequencing studies can be used to 
identify rare variants that may contribute to AF suscepti- 
bility. Recent sequencing studies screening lone AF 
patients for novel Cx40 coding variants have observed 
very different rare allele frequencies. In fifteen young, 
idiopathic AF patients of Western European descent, 
Gollob et al identified one germ line and three somatic 
nonsynonymous variants in the Cx40 coding region 
using DNA derived from atrial tissue samples and blood 
[9]. In contrast, we observed only one nonsynonymous 
variant in our 34 lone AF tissues. We hypothesize that 
this difference may be due to the high degree of selec- 
tion of lone AF subjects in the cohort studied by Gollob. 
In a prior analysis of Cx40 genetic variants using blood- 
derived DNA, Yang and colleagues reported the discov- 
ery of three nonsynonymous Cx40 variants in 218 lone 
AF patients of Chinese descent, and no novel variants in 
200 control subjects [26]. Our studies found no nonsy- 
nonymous Cx40 variants in blood-derived DNA samples 
from 91 lone AF subjects, which was not statistically dif- 
ferent from the frequency observed by Yang. Combined 
with our finding of no Cx40 nonsynonymous variants in 
30 mixed AF atrial tissue samples, these results suggest 
that Cx40 coding SNPs are uncommon in lone and sec- 
ondary AF patient populations. 

The current study has several limitations. First, all 
samples were predominantly derived from patients of 
European ancestry and the sample size was limited to 
91 lone AF blood-derived DNA samples and 67 atrial 
tissue samples, of which only 34 were lone AF subjects. 
Second, we did not possess paired blood and tissue 
samples, and thus were unable to determine whether 
the rare GJAS variants discovered in atrial tissue were 
of germline or somatic origin. We also did not have a 
large cohort of donor atrial tissue available that could 
be used to compare the frequency of rare variants in 
AF vs. control atria. Finally, our transfection study was 
performed in HL-1 cells, an atrial- derived, immorta- 
lized mouse cell line that demonstrates some atrial-like 
properties [12]. It is therefore possible that the Cx40 
3'UTR insertion/deletion polymorphism, which we 
observed to have no effect on reporter activity, could 
still modify gene expression in vivo in the human 
atrium. 
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Conclusions 

In 158 atria- and blood-derived DNA samples from sub- 
jects with AF, we observed only one nonsynonymous 
rare genetic variant in the GJAS gene. We determined 
that a 25 bp insertion-deletion polymorphism is com- 
mon in the GJAS 3' UTR, but we could detect no func- 
tional role for this variation. Our study suggests that 
aggregations of different rare nonsynonymous GJAS gen- 
etic variants are not commonly observed among AF 
patients. Although rare, GJAS variants that alter Cx40 
function may be directly causal of AF. 
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